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Abstract In critical percolation models, in a large cube there will typically 
be more than one cluster of comparable diameter. In 2D, the probability of 
>> 1 spanning clusters is of the order e~° ^ . In dimensions d > 6, when 
rj = the spanning clusters proliferate: for L oo the spanning probabil- 
ity tends to one, and there typically are ~ L^~^ spanning clusters of size 
comparable to |Cmaa;| ~ L^- The rigorous results confirm a generally ac- 
cepted picture for d > 6, but also correct some misconceptions concerning 
the uniqueness of the dominant cluster. We distinguish between two re- 
lated concepts: the Incipient Infinite Cluster, which is unique partly due to 
its construction, and the Incipient Spanning Clusters, which are not. The 
scaling limits of the ISC show interesting differences between low {d = 2) 
and high dimensions. In the latter case {d > 6?) we find indication that 
the double limit: infinite volume and zero lattice spacing, when properly 
defined would exhibit both percolation at the critical state and infinitely 
many infinite clusters. 
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1. Introduction 



l.a Main Results 



Some of the essential features of critical phenomena are reflected in the structure 
of the large correlated clusters observed in the critical state 0, ^ |] . For percolation, 
the relevant notion is of connectivity, and the characteristic feature of the criticality is 
the spontaneous formation of large connected clusters, of "macroscopic" diameter. In 
this paper we present rigorous results on the number of Incipient Spanning Clusters 
(ISC), in various dimensions. The term Spanning is applied generically to clusters which 
stretch across a large region, and connect different boundary segments (see Figure |I|, 
and the more detailed explanation in Section Incipient is the word we use for the 
tenuous structures seen at the critical models, at percolation threshold. 

Our goals are: i. clarify the issue of uniqueness, concerning which there has been 
some confusion, ii. present new bounds on the distribution of the number of ISC in 2D 
models, and iii. prove that in high dimensions Incipient Spanning Clusters behave in a 
rather different way than they do in 2D. 

More explicitly, we establish the following. 

1) In any dimension c? > 1, at p = pc-, the probability of there being more than one 
spanning cluster in a slab of the form [0, tL] x [0, LY~^ does not go to zero, as L ^ oo; 
the proof covers the case of t not too large, but independent of L (Theorem 2, below). 

2) For completeness, we add that above the critical point, p > pc, there typically is 
exactly one spanning cluster (for L >> 1) (Theorem Appendix C). Below the critical 
point there is none. 

3) In 2D critical models the number of spanning clusters is of finite mean, and its 
probability distribution satisfies: 



p I there are exactly k unconnected \>Ae"'^ fll) 
I spanning clusters in [0, L]^ I < e~"' 

where k counts clusters connecting the two opposite faces {xi = 0} and {xi = L} of 
[0, L]^, and the bounds hold uniformly in L (Theorem 3). 

The power k'^ seen in eq. ( p. . 1|) represents k'^li'^-^) , The analysis suggests a plausible 
argument for a more general validity of such lower bound in low dimensions, but it is 
not clear whether the actual rate is not slower for d=3,4,5. 

4) In dimensions c? > 6, assuming rj = [r] being the critical exponent seen in eq. ( p..4|) ): 



• the spanning probability tends to one (Theorem 4), and 
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• spanning clusters proliferate — 

Probp^ (the number of spanning clusters > o(l)L'^^^^ — > 1, (1-2) 

in the sense that the limiting statement applies to any substitution in which o(l) 
is replaced by a function of L which vanishes as L — oo. 

Furthermore, at least for clusters whose size (the number of points) is comparable 
to the maximal (ICImax), L'^~^ provides the actual rate of growth — not just a 
bound. 

• The size of the maximal cluster, as well as the size of the maximal spanning cluster, 
are typically 

l^maxl 5 |Cmax| ~ L , (1-3) 

in a sense made explicit in Theorem ^. 



The assumption on rj refers to the purported law: 

Const. 

T{x,y) = Probp^ [x sites x and y are connected) ~ -j " i^-'^J 



It is expected that 77 = above the upper-critical dimension = 6 |Q ||. Rigorous results 
in this direction were proven by Hara and Slade (reviewed in 0), who establish that a 
related, but somewhat weaker condition, holds at > 6 in models with sufficiently spread 
finite-range connections, and at somewhat higher dimensions for the nearest-neighbor 
percolation model. 

The results proven here for d > 6 reinforce the generally accepted picture, in which 
the proliferation of large clusters is related to the breakdown of hyperscaling (Coniglio 
||I|), and the "fractal" dimension of the large clusters stabilizes ai D = 4 (Aharony, 
Gefen and Kapitulnik 0, Alexander et.al. (AGNW) [[l^])- Our analysis is based on 
the diagrammatic bounds of Aizenman and Newman |]TT|. The high effectiveness of the 
method points to the validity of the perspective offered by AGNW |Ty] that for d > Q 
large clusters resemble randomly branched chains. 

The fact that the spanning probability tends to 1 has interesting implications for 
the nature of the continuum limit, mentioned in Section ^. 

Some of the topics discussed here are also the subject of the work of C. Borgs, 
J. Chayes, H. Kesten and J. Spencer [|12], where the properties of "incipient infinite 



clusters" are discussed within the context of an axiomatic description of the critical 
behavior. [The reader is advised, however, that the notions have not been coordinated. 
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and terms like IIC and ISC may be used in different ways. The terminology used in this 
paper is explained and motivated in greater detail in Section 5.] 

The prototype for the percolation models discussed here, and the systems to which 
we refer by default, are the nearest-neighbor bond, or site, percolation models on the 
regular lattice Z'^. However, with minimal clarifications the discussion applies also to 
percolation models with different short-scale characteristics, including finite range bond 
percolation, n.n. and next-nearest-neighbor site percolation, and also the continuum 
random dots models. To attain this robustness, the arguments employed refer mostly 
to large - scale connection events. The essential features (or at least those assumed 
throughout) are: independence for sufficiently separated regions [a fixed finite distance], 
and the vdB-K property which is presented in Section 0. 

Following are some less technical comments on the results. 



l.b Implications of the Multiplicity of Incipient Spanning Clusters 

The assertion that in 2D there is positive probability for more than one Incipient 
Spanning Cluster will not surprise mathematicians familiar with the theory developed 
by Russo [T^ and Seymour and Welsh [|14| . Nevertheless, even the 2D case of the general 



result (1) is in contradiction with statements which for quite a while have apparently 
been part of a wide consensus. That consensus was recently challenged, and corrected, 
by a report of a numerical work |T5| (see also |TB[)- As we attempt to make it plain 



(and as was already stated in ref. [[T^l)? the case of 2D is really simple. To prove the 
result for general dimension, we develop what may be regarded as a rigorous real - space 
renormalization group argument. 

The misconception of the uniqueness of the spanning clusters could have been 
facilitated by a number of factors: 

i. Scaling Theory 

The assumption that there typically is one dominant cluster (in a finite region [0, L^) 
helps in the explanation of the hyperscaling law, which is found to be obeyed in low 
dimensions. To re-emphasize that low level multiplicity (up to is consistent with 

hyperscaling, we present a version of the heuristic argument in Appendix A. It is seen 
there that the relevant condition is only that the number of large clusters grows slower 
than any power of L. 

ii. Uniqueness of the Infinite Cluster 

One of the results known under rather general assumptions (discreteness and "regular- 
ity") is that when there is percolation, the infinite cluster is unique |jl8|, [l^, |2^ (see 
however Sect. ^ below). The Uniqueness Theorems were occasionally quoted as an indi- 
cation that at the percolation threshold, where large-scale clusters are expected, there 
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is a single dominant cluster — referred to as the Incipient Infinite Cluster (IIC). Here, 
the convenient and otherwise appealing terminology might have suggested an incorrect 
deduction. 

iii. Terminology 

The term Incipient Infinite Cluster is too encompassing. The dominant cluster (s) in a 
large volume can be viewed from two perspectives, for which it may be better to have 
separate terminology. We propose to keep the term IIC for the large clusters as seen 
from the perspective of their sites (which presumably agrees with the conditioned process 
studied by Kesten in the context of 2D models). The same clusters, as viewed 
on the bulk (macroscopic) scale may be referred to as the Incipient Spanning Cluster (s). 
The distinction becomes more clear as one contrasts the microscopic and the macroscopic 
view of the system, and contemplates the relevant mathematical descriptions of the 
scaling limit — the limiting situation in which the ratio of the two scales diverges. A 
more detailed discussion of these issues, is presented in Section ^ which can be read 
in advance of the more technical Sections. Included in Section ^ are three possible 
definitions of the IIC, which still ought to be proven equivalent. For each of them it can 
be shown that if the percolation density vanishes at Pc then the IIC process exhibits a 
unique infinite cluster. The ISC in general are not unique (Section 2). 

Red — Blue bonds 

The multiplicity of ISC means also that the important notion of "red bonds" (Stanley 
^ 0, Coniglio |2^) needs some care. This notion is best explained in the situation in 



which a macroscopic piece of material is placed between two conducting plates which 
are at different electric potentials, and the bonds represent conducting elements. The 
"red bonds" are often interchangeably described as: i. bonds which are essential for the 
connection, i.e., are unavoidable for every path connecting the opposite boundary plates, 
or as ii. bonds which carry the full current. (Other elements in this classification are: 
"blue" - the backbone bonds, and "yellow" - the dangling ends.) As we now appreciate, 
these two definitions do not coincide. 

The non-uniqueness of the spanning cluster has the consequence that frequently 
the random configuration will not contain any bond which is essential for the connection 
(though there will be bonds essential for their specific spanning clusters). Interestingly, 
this does not mean that no bond carries the full current. If there are large differences in 
the total resistance of the distinct spanning clusters, than the bulk of the current may 
still be carried by just one of them. Its essential bonds meet the condition in the second 
definition of "red bonds" . However, as we see now, this terminology involves more than 
just the topology of the cluster. 

Despite the above complication, it is still desirable to have an elementary geometric 
concept expressing the fact that in typical configurations there are bonds whose removal 
(or blockage) will result in a change of the available connected routes which is visible 
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on a large scale R. For this purpose one may define a bond b to be essential on scale R 
( >> lattice scale) if there are two realized paths, of end-to-end distance R emanating 
from the opposite ends of the bond, which are not linked by any other path which stays 
within the distance R from b. 

Further discussion of the results is found in Section ^ 
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Figure 1: The geometric setups used in the discussion of the Incipient Spanning Clusters. 



Theorem 1 (Incipient Spanning Clusters occur) For any dimension d > 1, there is an 
L-independent (decreasing) function hd{t), which is strictly positive at least for < t < 
1/2, such that at the critical point for every < L < oo; 

/ the slab Sit = \0,tL] x \0,L]^'^-^'^ is traversed \ 
\ in the direction of the 1-st coordinate j 

Furthermore, hd{-) satisfies 

h{t) > 1 - Ae"™"^*- ast\0 . (2.2) 

The two dimensional case of this statement was derived in the works of Russo 
13| and Seymour and Welsh (conveniently summarized in ||2^), where it is also 



established that h2{t) > for all < t < oo. The result provides a very versatile tool for 
the study of 2D critical models (see below). We expect that also for arbitrary dimensions 
hd{t) > for all values of t, though that has not been established. The proof of Theorem 
1 is presented in Appendix B. It uses arguments which may be regarded as standard. 

Our main result for arbitrary dimension is the following. 
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Theorem 2 (There can be more than one ISC) For any dimension d > 1, there is a 
function Qdif), strictly positive for <t < to, with which for every finite size L: 

there is more than one 
DLit,p,)=Probp^\ I > gdit) . (2.3) 

spanning cluster in SL,t 



Remarks: 1) It should be appreciated that the bound in eq. ( ^.31 ) is satisfied only at 
p = Pc- Otherwise: 

lim DL{t,p [^Pc]) = (2.4) 

L^oo 

(for any t > 0). For p < Pc eq. ( p.4| ) is true because there are no spanning clusters (the 
connectivity function decays exponentially), and for p > pc it holds because there is only 
one spanning cluster (see Appendix C). 

2) As in Theorem 1, we expect the restriction to small t, not to be relevant for the strict 
positivity of gd{t)- I.e., if the probabilities of the corresponding macroscopic scale events 
are uniformly positive for some aspect ratio {to > 0) then that should also be true for 
any other value of t. However, we do not have a proof of such a principle for c? > 2. For 
2D, the strict positivity of g2(t) is directly implied by the RSW theory, as is explained 
in the next section. 

The proofs of Theorem 1 and Theorem 2 (for d > 3) have the flavor of a real - space 
renormalization argument. By considering suitably chosen local events, determined in 
blocks of scale L, it is shown that for any p: 



i. If RL,t is too small than the connectivity function decays exponentially, and hence 

P < Pc- 

ii. If [1 — RL/3,3t] + t is too small then there is percolation, and p > Pc- 

For p = Pc the term [1 — RL/3,3t] vanishes as t \ 0, and hence D^ t cannot be too 
small. 



The flrst argument is recapitulated in Appendix B. The second is given below. 

Proof of Theorem 2: As is explained in the next section, for 2D the RSW theory 



131 , p^] implies that eq. ( p.3|) holds with a function gd{t) which is positive for all t > 0. 

To prove Theorem 2 for ci > 3 we shall show that there is a constant b < oo (based 
on 2D considerations, and thus independent of d) such that for every p < Pc 

il-RL/3,3t)+DL,t > b (2.5) 
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for any L < oo. Since Ri^t > hd{t), that would establish the claim, with 

gdit) = b-[l-hd{3t)] (2.6) 

where limt^o[l ~ hd{t)] = 0, by Theorem 1. 

The following argument proves that if (for a suitable choice of the constant b) 
eq. ( |2.5| ) fails for some L than there is percolation in the slab [0, tL] x [0, L]'^^^ x R"^. Let 
us partition the slab into a planar collection of translates of the (i-dimensional rectangular 
cell [0,tL] X [0,L]'^~^ X [0,L/3]^. In a given configuration of the percolation model, a 
cell is declared regular if the following two conditions are met: 

i. the cell contains a spanning cluster (in the first direction), 

ii. the block made by joining 3x3 cells, with the given one at the center, is traversed 
by exactly one spanning cluster. 

Notice that if within the 2D array (of d-dimensional regions) there is chain of 
neighboring regular cells, then the spanning clusters of these cells belong to a common 
connected cluster. 

The last observation implies that there is percolation in the slab unless any finite 
region is encircled by a *-connected loop of cells which are not regular in our terminology. 
We estimate the probability of such events by a Peierls-type argument. 

While irregularity is not independent for adjacent cell, it is independent for cells 
with a gap of 2 in one of the two direction. Thus, it is convenient to estimate the 
probability of the existence of an encircling loop by focusing on sub-chains of minimal 
steps with such gap. There are 20 possibilities for each step in the chain. The Peierls 
estimate, suitably-modified, implies that a sufficient condition for percolation of the 
regular cells is 

Prob ( cell is irregular ) ■ 20 < 1 (2-7) 

We bound the probability of a cell to be irregular by adding the probabilities of 
the two possible failures. The result is that there is percolation unless eq. (|2.5|) holds, 
with b = 0.05. (That may not strike one as a very large number, but still it is positive 
uniformly in L). 

Since the probabilities Ri^t and D^ t are continuous in p (pertaining to regions of 
fixed finite size), if eq. (|2.5| ) fails at p then it does so also for p — e for e > small enough, 
and thus the conclusion is that either p > pc, or eq. (|2^ ) holds. 
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D 



□ 

□ 



Figure 2: The elementary mechanism for multiple spanning clusters in 2D critical models: 
a) strips traversed by connecting paths alternate with strips traversed by separating dual 
paths (which avoid the realized connections), b) long connected paths are formed through 
the intersection of elementary crossing events. 



We view the above argument as part of the rich legacy of the RSW theory. This 
term refers to a versatile method for construction which in combination with other 
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insights has been employed for a variety of far less elementary results: non-percolation 



at Pc 



26[| , support for the scaling theory 



sharpness of the phase transition 
(all the above for 2D models), some rigorous real space renormalization arguments 
and some extension to higher dimensions p8[. (Extension to the continuum 2D random 



Voronoy cell percolation model is worked out in [p9[|.) 

Our main goal in this section is to present large-deviation estimates for the prob- 
ability that the number of clusters spanning [0,L]^ exceeds (or equals) a large number 
n. 



3. b Bounds for the Number of ISC in 2D 

It is easy to be mislead about the large deviation asymptotics. By an application 
of the van den Berg - Kesten inequality 

KM < [RL,ir < e-li^'^^WI", (3.2) 

with /i2(l) = 1/2 in the self-dual case. Numerical results seem consistent with expo- 
nential decay in n (of the form e~'-^°"''*- ") Nevertheless the decay rate is different 
(faster). 



Remark: The inequality of van den Berg and Kesten states that for independent 
percolation (or, generally, systems of independent variables) the probability of the dis- 
joint occurrence of two, or more, events (e.g., the existence of two separate connecting 
paths) is dominated by the product of their probabilities. The statement is not obvious 
since disjoint occurrence does not refer to a-priori specified separate regions. 

The faster-than-exponential decay is caused by the mutual exclusion, which limits 
the clusters to reduced regions. (For the lower bound seen in eq. (|3.1| ) the clusters are 
produced in disjoint strips). 

Theorem 3 In planar percolation models, at p = pc the probability that [0,L]^ is tra- 
versed ("left ^ right") by at least n separate clusters satisfies [for n » 1] 

K,{n) - , ^ , (3.3) 

I < e"'^ " 

where the bounds involve different positive constants. 

Remark: Prob{ exactly n crossings ) behaves just as Prob{ at least n crossings ) = 
Klin), since the former is KL{n) — KL{n+l), and (by the vdB-K inequality) KL{n + l) < 
KL{l)-KL{n) {=KL{n)/2 ) . 
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Proof: The lower bound in eq. ( p.3| ) is obtained through the construction discussed 
above, supplementing eq. with the known fact that h2{t) > [^]e^°'^ (proven by 

constructing a chain of 1 x 2 brick crossing events; as seen in Figure ^.b |jl3|, [1^, ^ [28[1). 

The upper bound is slightly more involved. Let us first ask how many distinct 
spanning clusters will be seen in a tall strip Lh x {Ly > Lh)- It is natural to expect 
that the number, call it A^l^.l^ will typically be of the order of Ly/Lh, since, roughly, 
there is a finite number of distinct crossings in each of the squares comprising the tall 
stack. This motivates the following bound, which is proven below using Lemma |l]. 



Prob. 



Pc 



< e 



-f{u) L^/Lh 



(3.4) 



with f{u) > for u large enough. 

For a moment, let us assume eq. (|3.4|) . For our analysis, we fix m at a value for 
which f[u) > 0; the optimal choice being one maximizing f{u)/u'^ (one may note that 
f{u) = for M < 1/2). 

To estimate the multiple-spanning probability KL{n) (for u < n < L),we partition 
the L X L square into n/u vertical strips, so that the height /base ratio satisfies n = 
uL/Lh- If the square is spanned by at least n disconnected clusters, then that is also 
the case for each of the vertical strips. (Though the converse need not be true.) These 
being independent events, the probability is bounded by: 



KL{n) < 



n/u 



/{") 2 

< e 



which concludes the proof of eq. (p.3|) . 

To prove the large - deviations estimate eq. (|3? 
related statement. 



(3.5) 



we first derive the following 



Lemma 1 



(3.6) 



with (j){t) < oo at least for t small enough. 

[A - posteriori, we'll find (j){t) to be finite for all t] E( ) represents expectation value. 



Proof: To avoid more complicated analysis, we shall take advantage of a convenient 
shortcut. For each length L, let Nl be the (random) number of distinct connected 
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clusters connecting the left face of a square with the vertical line which includes the 
right face, allowing the paths to meander into the entire vertical strip containing the 
square. We claim that there is a function 0(t), finite for t small enough, such that: 

/■oo _ 

Ep.(e* = 1 + t e*" ProhpXNi > x) dx 

^ 

< e'^(*) . (3.7) 

The claim is based on two observations: i. Proh{NL >k)< Prob^Ni > 1)^ (using 
the vdB-K inequality), and ii. Prob{Ni > 1) is uniformly < 1 (a simple construction 
using the RSW theory). 

A moment of reflection shows that the random variable Nr^ r is dominated, in the 
stochastic sense, by the sum of L^/Lh independent copies of N^^ [that is the short-cut.] 
To see this, let us partition the strip's left vertical edge (of height L„) into intervals 
of length L/i, and for each integer 1 < k < L^/L^ let be the number of spanning 
clusters whose lowest intersection with the strip's left edge falls within the fc-th interval 
(counting from below). 

The conditional distribution of Nk, conditioned on {Nn : n < k }, can be studied 
by further conditioning on the exact location of the previously counted clusters. The 
net effect of the conditioning is an excluded volume. Regardless of it, the conditional 
distribution of is dominated by that of N^^, since the latter is recovered through the 
addition of two connected regions (one above the free region and one below) in which 
further clusters may be found. We rely here on the fact that in 2D, the extra connections 
cannot diminish the number of spanning clusters by joining different ones, since they 
can reach only the "lowest" (or the "highest" ) cluster. The most relevant implication of 
the stochastic-domination relation is: 

as claimed in Lemma |I|. (I thank H. Kesten for correcting the first draft's definition of 

The remaining path from eq. (|3.8|) to eq. ( ^^ involves the Chebyshev estimate: 
Prohp^ (Nl„l. >u^) < Ep,(e* ^^'^■^") e"* " ^"/^'^ < , (3.9) 

which shows that eq. (|3.4| ) holds with f{u) the Legendre transform: 

f{u) = sup[tn - 0(t)] . (3.10) 
t 
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(3.8) 



Figure 3: Two multiple-crossing events which are related by the map z j-gSTr^/L ^^j^]-^ 
R/r = e^'^^ and L, r >> lattice spacing). Qualitatively, our bounds are invariant 
under such conformal maps. 

Before we turn to the discussion of the spanning clusters in high dimensions, let 
us note that the arguments used to derive Theorem 3, provide also information on two 
other questions (depicted in Figure |). The first relates to the rate of exponential decay 
for the probability that a long strip ([0, fcL] x [0,L]) is traversed by n (< L) distinct 
spanning clusters. We find that the rate of decay in the aspect ratio k (the "mass" in 
field-theoretic jargon) is of the order of n^: 



where the focus is on >> 1. The proof is just the argument of Theorem 3, carried for 
rectangles rather than squares. 

The second statement concerns the power law for the decay of the probability that 
an annulus is spanned by n distinct clusters. This time, it is the characteristic exponent 
which depends on n: 




the strip [0, kL\ x [0, L] is traversed (in 
its long direction) by n distinct clusters 




(3.11) 
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Prob ( annulus {xG-R :r<|X|<_R}\i ,^ 

\ is traversed by n distinct clusters / ' - \ ■ J 



> A (r/i?)-""V(2-) 



where it understood that 1 < n < r/ lattice spacing. As explained in Figure 0, the events 
seen in eq. and eq. ( p.l2|) are related by a conformal map. This map indicates also 



how to transcribe our proof of Theorem 3 into a proof of eq. ( p. 121 ) for an annulus with 
a cut. It is also easy to adapt the argument to the full annulus and, correspondingly, to 
the rectangle with periodic boundary conditions. 

The constants a with which the above equations are proven are not identical, 
though it is natural to conjecture that the asymptotic values (in the sense seen below) 
are the same, and agree with the limit: 

lim lim ^logKiin) (= a(^asymp.) ?) • (3.13) 

It seems plausible that the exact value could be determined by methods based on the 
Coulomb gas representation and / or conformal field theory [^|. Further discussion 



of the purported conformal invariance of the scaling limit can be found in references 



4. Above the Upper Critical Dimension 

J^.a The Relevant Condition 

We shall now prove that in contrast with the situation in 2D, above the upper 
critical dimension spanning clusters are sure to occur, and they proliferate in the way 
described in the introduction. A convenient sufficiency condition is that: 

d>6 (4.1) 

and the connectivity function defined by eq. ( |1.4|) (at p = pc) satisfies: 

T{x,y) < Const./ |x-y|('^-2+r,) ^4^3) 

with 

r/ = 0. (4.3) 
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In this article, the relation (<) seen in eq. ([4 .21) means that 




with a pair of possibly distinct constants, < C" < C < oo. 

It is generally expected that d > 6 =^ rj = 0. As mentioned in the introduction, 
there are mathematical results offering partial support to this claim, however some gaps 
remain |^. The statement which is established in ref. |0| for high dimensional models 
is slightly weaker than eq. ( [4.2|) and eq. (|4.3| ), but sufficient for the "triangle condition" 
of ref. [0 ( " V-diagram" (pc) < oo). A strategy which was successful in the derivation 
of some other results about the critical behavior in high dimensions was to base the 
analysis on such a weaker assumption (proving "V-condition" =^ 7 = 1 |TT|, (3 = 1 and 
6 = 2 We do not pursue that track here. 



4-b The Boundary Conditions 

The number and the sizes of the connected clusters in a finite volume are affected 
by the boundary conditions (b.c). (I thank C. Borgs and J. Chayes for calling my 
attention to this point.) The possibilities include: 

• Free b.c. — sites can be connected only by paths within the specified volume, A. 

• Bulk b.c. — the given region A is viewed as part of a much larger domain, taken 
to be all of R'^. Sites in A can be connected by paths which stray outside of that 
region. 

• Periodic b.c. — in the well familiar sense. 

(Omitted here are the Wired b.c, for which Nwired = !•) 

The obvious relations for the numbers of [left-right spanning] clusters, are: 

Npree > Neulk , Np^r (4.5) 

with opposite inequalities for the maximal cluster size. 

The terminology we adapt here is based on the Bulkh.c. This choice facilitates the 
presentation of the main ideas, given that our starting point is an explicit assumption 
concerning the two point function in the bulk. With additional information on the two 
point function, the method should be applicable also to the other b.c. 

Question: In a sense made more precise below, we shall show that the maximal cluster 
size (Bulk b.c.) scales as |C|max ~ L"^- Concerning the other boundary conditions: we 
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expect that the Free b.c. scahng is similar, in hne with the prevaihng theory 0], but 
we venture the guess that the Periodic case is significantly different, possibly scaling as 
|C*|max ~ iJ''^!'^ . The latter is the scaling law for the complete graph [38| (see also |12]). 



Our guess is motivated by an analogy with a phenomenon which has been noticed for 



the (p2 field theory above its upper critical dimension Is the guess correct? It will 
be of interest to see numerical work on this question. 

J^.c The Number (L'^^^) and Size (L^) of the Spanning Clusters 

The geometric setup we use is that of the d-dimensional cube, denoting A = A^, = 
[~L,LY, and <9A_(+) = {x e A^ | xi = -L(+L)} (the "left" and "right" boundaries). 
The symbol C would range over the connected clusters in A (connected in the "bulk 
b.c." sense). 

Spanning clusters are defined now as those intersecting both 9A_ and 9A+. It can, 
however, be easily seen that none of the results will change (except for constants) if that 
notion is modified by declaring as spanning only those clusters which reach each of the 
2d boundary faces. 

The number of spanning clusters (which is a random variable) is denoted by N^. 
We shall also be interested in the number NLy^l of clusters spanning A^^ which are visible 
in a much smaller window Bw = [—W^ with 1 (=lattice spacing) << W « L: 

Nl,w = card {C C A^ | C n 9A_ ^ 0, CD dA+ ^ and C n 5vi/ ^ 0} • (4.6) 
(the location of the window within A will not play an important role.) 
The main results in this section are the two Theorems found below. 



Theorem 4 In any dimension d > Q, if eq. ([j.Sj holds with rj = 0, then, as L oo: 

1 ) The spanning probability (at p = pc) tends to 1. 

2) The probability of observing a spanning cluster within the window = [—W,W]'^ 
also tends to 1, provided 

W 

^ oo . (4.7) 



^2/(d-4) 

L— >oo 



3) The numbers of spanning clusters (in A and those observed in By/) satisfy: 

Probp^ [Nl > o(l)L'^-^) 1 . (4.^ 



L— >oo 
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Probp^ INl^w > 1 • (4.9) 

L, W{<L)^oo 

in the sense that the statements are true for any function of L (and W), denoted here 
0(1), which tends to as L,W ^ 00. 



The proliferation of the spanning clusters bears some relation to their low ("frac- 
tal") dimension. It is often stated that their dimension reaches D = 4 and remains 
at that level, as the lattice dimension is increased over the upper critical value d = Q 
(1^, [1^, |I|). The following statement offers some rigorous support to this claim. 

Theorem 5 In d > 6 dimensions, assuming the power-low behavior (eq. U-^ )) with 
r] = 0, 

(r/ M < c log L ■ . 

max[(^P)l \C\ ~ — ^ 1 (4.10) 
'^^^ > 0(1) ' 



L—*oo 



and 



Probp^ 



V 



r/ M < c logL ■ , 

max \C\ndAL ~ — ^ 1 (4.11) 
> 0(1) 



L— ►oo 

where [(sp)] represents an optional restrictions to spanning clusters, and o(l) is to be 
interpreted as in Theorem^. 

In the derivation of these results we make use of the diagrammatic bounds of ref. 
TT| . The most elementary of these is the tree-diagram bound on the k-point connectivity 



function, which directly yields: 

Lemma 2 In any dimension, assuming the power-low behavior eq. I \4-^ ), for any k >2 
E,,(|C|L.) < E,J ^ \cA < fc! C,^- ■ (4^12) 

\CcAl / 
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and 



Ep. I I] |Cn5iy| |Cn(9A_| |Cn9A+| I < Const. W'^L'^-^'^ . (4.14) 



Somewhat less immediate are the lower bounds seen in the following auxiliary 
statement. 



Lemma 3 In any dimension d > 6, if eq. l \4-3i ) holds with rj = 0, then for k = 1, 2 

Probpi \C n dA^\ \C n dA+\ ~ ^^"^^^ x L'^'^ L^'^ \ — > 1 (4.15) 
IccAl > o(l) 



L— >oo 



and 



Probpl \C n Bw\ \C n dA-l'' \C n dA+\ ~ ^^"^^^ X VV^lA — ^ 1 (4.16) 
\cc\l > o(l) / 

L^oo 

[where we read W'^L'^ as * L^]. 

Before deriving the two Lemmas (in the next subsection), let us go over the de- 
ductions of the Theorems (| and ^) , from Lemma ^ and Lemma |^ . 

Proof Theorem ^ {assuming Lemma ^ and 

We deduce that |C|max = niaxc^Ai IC" H A^l is typically not much larger than L"^, from 



the rate of growth of the moments described in eq. ( [4.12| ) and eq. ( [4.13| ), combined with 
the Chebyshev inequality (see eq. ( [4.30| ) below): 

E( \c \^ 

Prob,^ (|C_| > a{L) • L^) < ^^^ym^ ■ (4.17) 

A comparison with eq. ( [4.12| ) shows that a natural choice is A = (4 — 2t]) and a{L) = 
alogL. Optimizing over k (chosen so that k ^ a{L)/Cd ), one learns that 

ProV {\Cma.\ < alogL- L^^'^")) 1 , (4.18) 

L^oo 
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provided a > Cd{d — Q + 3ri). This proves the upper bound on \C\max claimed in [4.10 
in a somewhat more general form. A similar use of the Chebyshev inequality leads from 
eq. ( 4.13|) to the conclusion that typically maxc^A^ |C fl (9Al| < a' log LL^~^^, as claimed 
in eq. (^Tl\ ). 

To deduce that \C\max is typically not much smaller than the above upper bound 
(L^), we use : 



Id > \C(^P)\ > |LIIC>A_| \HWA+\ 



E|CnA| 


\C n 9A_ 




CndA+\ 


E\Cr\dA^\ 


\c n dA+ 





(where the sums are over the connected clusters of A = Aj^) combined with the infor- 
mation about the typical values, provided by eq. ( [4.15| ) with k = 1, and eq. ( [4.1(j| ) at 
W = L. Likewise, the lower bound claimed in eq. ( 4.11| ) on max^^A^ '^^^^l |C| H dA^ 
follows by comparing the lower bound in eq. ( [4.15| ) for k = 2 with the upper bound in 
the same equation for k = 1. 

Proof Theorem ^ (assuming Lemma H and To show that typically there are at least 
o{l)L'^~^ spanning clusters we use: 



E\Cn dA_ 




CndA+\ 


max 


\c n dA. 




\CndA+\ 



Nl > ' 1^ — ^. ' ,^ — — ^ . (4.20) 



Typically, the numerator is > o{l) L'^ (by Lemma ^ while the denominator is 
< l/o(l) L^-^ (by Theorem 1). Hence Nl > o(l) L'^ ^, as claimed in eq. 



At the risk of sounding too repetitive, for the number of the spanning clusters seen 
in the window Bw we use: 



E\CnBw 




CndA_ 




CndA+ 




max 


cnBwl 


max 


CndA 






CndA+\ 



which leads to the conclusion that typically: 

y/d^A yyd-4 

Nl,w > o{l)^^, = 0(1)^ . (4.22) 

Under any conditions which guarantee N^^^y > 1, we are assured that a spanning cluster 
intersects the box B^. This concludes the proof of Theorem ^ 

A particular conclusion which seems worth emphasizing is: 

, j the window [—W, W]^ is connected to the boundary 1 , , 



provided a < {d — A)/2 (assuming all along that d > 6, and t] = 0). 
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4-d Derivation of the Moment Bounds 

We start from the elementary but important identity 

e(E = ^(y. i[xieC] ...iKgC]) 

\ C Xi_,...,Xk / 

= E E I[ xi, . . . , Xfc belong to a common [spanning] cluster] j 
= r{xi,---,Xk) (4.24) 



where ![-] is an indicator function and Tfc(xi, ...,xa:) is the probability that the k points 
are all connected. 



B 



{in} - 



T(A,B) T(C,D) 



C D 




Figure 4: The inequality used for estimating the variance of I] |Cp. The trunca- 
tion adds in the diagram two vertices and three lines. The result is an extra factor 
of Const. /,2dy'^3(d-2+r7) _ Qg^gf^ / L'^^^^^'^ , which abovc the upper critical dimension 
tends to for large separations. 

Next, the argument will employ diagrammatic bounds on the connectivity func- 



tions, following the approach presented in ref. |Tl|]. The technique has been further sim- 
plified through the development of the van den Berg - Kesten inequality |3l|] mentioned 
in Section 3. There are two principles whose repeated application leads to Lemma ^ and 
Lemma ^ (We shall attempt to avoid obvious repetition of similar arguments.) They 
are: 
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1) the tree-diagram bound |TT|, and 

2) the truncation lemma (Lemma proven below. 

The tree-diagram bound |Tl| states that Tk{xi, . . . ,Xk) is dominated by the sum 
of products of the two point function t{u, v), arranged in the form of tree diagrams with 
the external vertices xi,...,Xk- Its intuitive explanation is that in order for the given 
k sites to be connected there has to be a connecting tree. The vdB-K inequality |^ 



(or the alternative argument which was originally used in ref. fTTj]) permits to bound 



the probability by the sum over the corresponding tree diagrams. (The bound is not 
totally obvious since it involves sum over tree diagrams, and not tree graphs, i.e., all 
trees embedded in the lattice. The latter sum is much too large.) 

For k = 3 the inequality is: 

r3(xi,a;2,a;3) < ^ r(xi, m) r(x2, m) r(a;3, m) . (4.25) 

u 

For general k > 3 there are {2k — 5)!! < 2^A;! tree diagrams for each set of the external 
vertices, and each diagram has {2k — 2) vertices (external + internal) and {2k — 3) lines. 
For example, for the quantity seen in eq. ( [4.24| ) the sum yields: 



A different argument is needed for the lower bounds seen in Lemma ^ These rely 
on a more delicate estimate of the difference between two comparable terms, which under 
the right conditions cancel up to a small remainder. Following is that auxiliary result. 



Lemma 4 For independent percolation on any graph, 
I 



< Prob 



Tr 



,{Xi, ...,Xn)x Tn{yi, ■ ■ ■ .Vn) < 



{xi, . . . , Xn] are all connected, and 
{yi, . . . , y^} are all connected 
< Proh { xi, . . . , Xn, and yi, . . . ,yn belong to a common cluster) , (4.27) 



V 



and, for n = 2: 

I 



< Prob 



t{xi,X2) XT{yi,y2) < (4.28) 



{xi,X2} and {2/1,2/2} 
are pairwise connected 
— '^'''i^i^'^) ''"(z/i)^) t{v,u) t{v,X2) r(t>, ^2) + [yi ^2/2 permutation 



23 



Proof: The positivity of the difference seen in eq. ( ^.27| ) is obtained by the standard 
monotonicity argument (the Harris - FKG inequahty 0): the probabihty of the 
simultaneous connection of the two pairs is greater that the product of probabihties. 

In the opposite direction, if the two collections of sites are interconnected, then 
they are either in disjoint clusters or in a common cluster. The probability of the first 
alternative is smaller than the product which is subtracted (by the vdB-K inequality 

Thus, we are led to eq. ( |4.27D . 

We could in fact make a stronger statement: by separating not just the case that 
the two sets are in disjoint clusters, but the case that they are disjointly connected, we 
are left with a smaller remainder. In each of the configurations contributing to it there 
is a tree subgraph connecting all the 2k vertices. The tree is bound to have a link (a 
pair of internal vertices {u,v}) whose removal will split it into two subgraphs of equal 
numbers of vertices, and there is a further constraint that this link cannot separate the 
x— sites from the y-sites. For n = 2, that leads to eq. (|4.28|). 



Proof of Lemma g (following ref. |[TT[| ): 

The first claim, eq. ( [4. 121) , is proven in eq. ( [4.26|) . The other statements of Lemma ^ are 



derived similarly; e.g., for eq. ( |4.14| ) one may start with: 



Epe ( E I'^^^vfI |Cn(9A_| |Cn<9A+| ) = J2 



Ts{x,y,z), (4.29) 

xeBw, y&dA-, zedA- 



and follow by an application of the tree-diagram bound, eq. ([4.25|) . 



Proof of Lemma ^: 

Let us consider first eq. ( [4.15|) . The upper bound seen there is a direct consequence of 
the value of the mean and the Chebyshev inequality: 

Prob{X >t) < . (4.30) 

A somewhat more delicate argument is needed to deduce that typically the sum is not 
much below the mean-value level. 



As in |Cp version of the identity eq. ( ^4.24|) , it is easy to see that: 



^ |C n dA-\ \C n dA+\ = l[x andy are connected ] := K , (4.31) 

CcAl xe9A_, yedA+ 

where K is defined by the sum. The mean value of K is easily determined: 

xg9A_, J/G9A+ 

< Const. L^id-^^ / L'^-^+v = Const. L"^-". (4.32) 
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L 



Figure 5: The dominant tree-diagram in the bound on Var(i^). As in Fig. 3, the 
truncation adds to the graph two extra sites and three extra hues, resulting here in the 
multiphcative factor Const. /W'^~^~^^'^ . 



The corresponding expressions for the second moment, and for the variance, are: 
E(K^) = ^ Prob (xi is connected to yi, and X2 is connected to y2 ) 



and 

Var{K) 



Xi,X2 G 

yi,y2 e A4 



Prob 

xi, X2 e A_ 

2/1,2/2 G A+ 



{xi,yi} and {0:2, 2/2} 
are pairwise connected 



(4.33) 

(4.34) 

T{xi,yi) X r(x2,2/2) 



Using eq. (|4.28|) of the truncation Lemma, Lemma ^ we obtain the following 
estimate (where E should be read as EpJ: 



E 



K 
E{K) 



E{K^) - \E{K)\' 
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< Const. L^VL^^^-^") = . (4.35) 

In essence, the "truncation" results diagrammatically in the addition of two sites and 
three lines, which translates to the multiplicative factor of order L^'^/L^*^'^"^"*"'') = 
^d-6+3r) Under the right conditions, eq. ( |4.35| ) implies that only very seldom will K 
differ by a significant factor from the mean E{K). 

We can proceed only under the assumption that d — 6 + St] > 0, for otherwise 
the last bound still allows the typical values of KfElK) to be arbitrarily small. That 
however is ruled out if the right side of eq. ( [4.35|) is o(l), in which case there is a constant 
b > (determined through eq. ([4.32|) ) for which 

Probp^ (k > b L^^+^-"^^ — y 1 . (4.36) 



Since J2^^^^ |Cp > K that implies the lower bound claimed in eq. (|4.15|) , and thus 
concludes its proof. 

The derivation of eq. (|4.15| ) is very similar, with K replaced by K defined by the 
following modification of eq. ( [4. 311 ) 

|Cni?vK| |Cn9A_| |Cnc?A+| = ^ I[a;, and 2; are connected ] := K. 

^^^^ xeBw 

y e 9A_ 
z e dA+ 

(4.37) 

In the study of the variance Var{K), the role of equation eq. (|4.28|) is taken by eq. ( [4.2?| ) 



combined with the tree diagram — bound. As before the "truncation" results diagram- 
matically in the addition of two sites and three lines (see Figure^), but now the main con- 
tribution comes from diagrams with two internal sites within distance of order W from 
the window B^. That translates to the multiplicative factor Const. W'^'^/W^^'^~'^^'^^ = 
Const./ 

Remarks: 1. Purists may note that the above proof, that spanning is certain, actually 
requires only the weaker condition: 

d-6 + 3ri>0 , (4.38) 

under which the conclusion is that the number of spanning clusters is typically greater 
than o(l)L'^~^+^'^. This can be turned around to say that in dimensions in which the 
spanning probability does not tend to 1: rj < The importance of this improvement 
is dimmed by the fact that according to numerical estimates 77 < for c? = 3, 4, 5. 

2. It is natural at this point to enquire about an upper bound on the number 
of spanning clusters. Comparing eq. ([4.12|) with Theorem ^ we see that for counting 
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clusters whose volume is comparable to the maximal the lower bound we got is sharp up 
to xL°^^\ However our estimates will not detect the presence of a large number (larger 
that L'^"^, but not too large) of sufficiently thin spanning clusters. 



5. The lie, the ISC, and the ScaUng Limit 

5. a The Microscopic and the Macroscopic Perspectives 

In this section we place the above results in the context of the interplay between 
the different scales on which the system can be viewed. 

While percolation models are often presented on an infinite lattice, many of the 
interesting questions relate to finite, though large, systems. There are therefore at least 
two scales of interest (we focus on two but it will also be interesting to bring in a 
third, intermediate scale, which could allow a better mathematical expression of the 
renormalization group approach): 

i. The microscopic scale, for which the unit length (a) is the lattice spacing, or the 
size of the dots in the random dot model. This is the scale on which the elementary 
connections are defined. 

ii. The macroscopic scale, for which the unit length could be the system's width (L). 
This is the scale of the clusters on which we focused. Of particular interest are the 
clusters which connect two different boundary faces. We generically refer to those 
as spanning clusters (in different geometries, see Figure 0). 

The situation in which L » a appears rather differently from the two perspec- 
tives. 

On the microscopic scale the entire system appears vast, and it is mathematically 
advantageous to consider it infinite, and homogeneous (in the suitable sense). The 
relevant limit is L ^ oo. The local structure of the very large clusters can be studied 
on that scale but only if one first centers the viewing field on one of their sites. That 
requires some effort, since at p = pc the density of the union of macroscopic-scale 
clusters vanishes. Thus is born the Incipient Infinite Cluster, which can be defined quite 
precisely through a number of procedures, most of them presumably yielding the same 
limiting object (see below). 
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On the macroscopic scale, the situation Lja » 1 is expressed through the scaling 
limit a ^ 0. This hmit lends itself naturally to the discussion of the higher symmetry, 
including possibly conformal invariance, which seems to emerge on the large scales in 
the critical regime. 

When we count the number of the spanning clusters, or other clusters of linear 
extent comparable with the system's size (say diameter exceeding 0.8L, where 0.8 can 
be replaced by any fixed < s < 1), we are obviously taking the macroscopic perspective. 

One should avoid being too dogmatic, as some of the features of the large clusters 
appear similar when one looks either down from the macroscopic scale, or up from the 
microscopic scale (e.g., the structure of voids, discussed in the recent study ref. p2| ). 



However, ignoring the above scale distinctions one runs the peril of pitfalls and para- 
doxes. 

An illuminating example is the following delightful paradox due to D. Stauffer 
[^. There is an early rigorous result [|r^ (even preceding the proofs of the uniqueness 
of the infinite cluster [jl9|, ^) that in an infinite translation-invariant system (meeting 
a certain regularity condition) the number of infinite clusters is almost surely constant: 
either 0, 1 or oo. In a finite volume, the term "infinite clusters" ought to cover the 
spanning clusters. Since infinity is ruled out, and we know that at p = pc spanning 
clusters do occur, we are led to the conclusion that typical configurations of critical 
models, exhibit a unique spanning cluster. That however is in direct contradiction with 
Theorem 0. 

The answer to Stauffer's paradox lies in a careful analysis of the relevant meaning 
of infinity. Positive probability of there being an infinite cluster (on the short scale) 
means that for some fixed K (large enough) the probability for observing a cluster of 
size R within the distance K from the lattice origin (or alternatively from an "averaged 
site" in the large box) does not tend to zero in the limit R,L — »• oo. (For independent 
percolation one may take R = L; the distinction is needed for correlated percolation 
models where "thermalization" may require the limit L — oo to precede i? — >■ oo.) We 
leave it as an exercise to the reader to determine the finite-volume translation of the 
statement that there are more than two infinite clusters in the infinite volume limit. 

At p = Pc no infinite clusters are seen in the infinite translation-invariant limit, 
assuming that the percolation probability vanishes (as is proven for 2D and for suffi- 
ciently high dimensions). Spanning clusters and other macroscopic-size clusters occur, 
but the density of their union tends to zero, and thus they leave no trace in the usual 
infinite volume limit. (Their microscopic scale separation also increases, and the IIC 
construction yields a unique infinite cluster.) As we see, the paradoxical contradiction 
might be first "corrected" as not between one [infinite cluster] and a number greater 
than that [of spanning clusters], but between zero and a number greater than 1. The 
resolution is through a finer appreciation of the difference in the scales on which the 
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objects are observed, and defined. 

That leaves us with another paradoxical observation: the infinite lattice represents 
just a small part of the finite system. This fact has been noticed before [Q, when it 
was appreciated that for ferromagnetic spin systems with the order parameter M = 
lim^j^oo I < o'oax > \^^'^ a distinction ought to be made between the Short - Long - 
Range Order parameter, defined on the infinite system, and the Long - Long - Range 
Order parameter for which the distance between the spins scales along with the box 



size. (Ipso facto, for the case considered in the two order parameters can nowadays 



be proven equal, but the distinction is of relevance for expectation values of quantities 
which are somewhat less local than the spins, as the Stauffer paradox demonstrates.) 

In mathematical terms, a careful inspection of the notion of measurability in prod- 
uct spaces also shows that the standard infinite system (e.g., the lattice) is in essence 
just the collection of all regions which remain at fixed microscopic-scale distances from 
the point represented by the lattice origin. Thus, again: the infinite system represents 
just a small part of the finite region studied in simulations. 

Let us add that it is an interesting question of what mathematical framework can 
be used to contain the macroscopic information of relevance in the limit a — > 0. The 
spanning clusters and other macroscopic-scale clusters have the appearance of fractal 
objects. However, it turns out that the standard description of fractals does not capture 
the relevant information on the realized connections. Thus a new formalism is required. 
A proposal for the limiting description of the percolation "Web" is presented in and 
work in progress. It is also of interest to take the double limit: a — > and L —>■ oo, 
which is needed if one wants to ask about percolation in the scaling limit, on which more 
is said below. 



5. b Conjectured Properties of the IIC and the ISC 

Of the different constructions which were proposed for what one may call the 
Incipient Infinite Cluster, our notion is closest to that of Kesten I^T] (another proposal 
P3[ is mentioned below). The goal is to provide the mathematical description of the 
large clusters, down at the microscopic scale. 

One may guess that the following three algorithms for the construction of a prob- 
ability measure of the percolation model in which a chosen site is constrained to be 
connected to infinity, will have common limits. The infinite cluster seen there would be 
called the IIC, and the limiting measure will provide its probabilistic description. 

a. Generate the probability distribution for a random environment by a two step 
procedure. First pick a typical random configuration, and then center it relative 
to one of the sites on the large clusters, sampling with equal weights over all the 
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sites connected to the boundary of [0, L]'^, possibly with a corrective exclusion of 
the sites closer to the boundary than some g{L), with g{L) — > oo for L ^ oo (e.g., 
g{L) = O.IL). 

b. Take the conditional distribution conditioned on the origin being connected dis- 
tance L away, and let L oo. 

c. Raise p over pc, condition on the origin being in the infinite cluster (which now 
exists), and let p\ Pc- 



For 2D, Kesten proved the convergence and equality of the limits b. and c. 
Presumably, the method should cover also a. (which seems to be the most efficient for 
simulations). While the result was not yet extended to higher dimensions, it is expected 
to be true for general d > 1. Regardless of that, under the assumption Poo(Pc) = one 
can prove that in each of the constructions, the limiting measure will exhibit exactly 
one infinite cluster (for p = pc). This uniqueness, however, is not the consequence of 
any of the uniqueness results cited above, since the measure (s) are neither regular nor 
translation invariant. 



Another construction of IIC was suggested in ref. where an infinite cluster 
is formed by enhancing the the bond densities at a rate which tapers off as ||x|| —>■ 
oo). Presumably, the different constructions will agree in the further limit a —>■ (the 
continuum limit). 

One may be puzzled by the discrepancy between the multiplicity of the Incipient 
Spanning Clusters (non-unique by Theorem ^ and the uniqueness claimed for the IIC. 
The long answer was provided above. The short answer is that the different ISC are at 
increasing distances apart (on the lattice scale). 

We have also a natural conjecture concerning the scaling limits of the Spanning 
Clusters (not just at Pc). We guess that for any < p < 1 one would find just one of the 
three possibilities, which would be equally valid, with probability which tends to one, 
for any bulk shape of the type seen in Figure ^ 

i. No clusters attain size visible in the scaling limit. 

ii. There is a unique macroscopic connected cluster filling the region densely (in the 
macroscopic-scale sense). 

iii. Multiple clusters with macroscopic-scale diameters will be seen in any rectangular 
region (regardless of its aspect ratio). 

Existing results permit one to easily deduce: 
p < Pc =^ alternative i. 
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p > Pc =^ alternative ii. (see Appendix C) (5.1) 
for 2D: p = Pc =^ alternative iii. (see Section 3) 

and in Section 2 it was established that 

any d > 1: p = Pc =^ a restricted alternative iii. , (5.2) 

where the restriction is to rectangular regions with a conveniently small aspect ratio. 
Since the discussion refers to the scaling limit, it seems most natural to expect the 
restriction to be irrelevant in any dimension, though that is not proven here. Such a 



result will amount to a significant extension of the theory of Russo and Seymour 



and Welsh 14 , which is based on intrinsically 2D arguments. 



5. c Distinction between Type I and Type II Models 

We shall not discuss here the interesting question of what mathematical object will 
provide a suitable description of the structure emerging in the continuum limit. Let us, 
however, point out that a distinction ought to be made between two classes of critical 
models. We characterize them as follows. 



Type I models: The function 

1 f the set \—sL,sLY is connected \ r , ^ „x 
hm sup ProV boundary of [-L, L]^ ) = ^^'^ ^^'^^ 

is strictly less than one, for some < s < 1; in which case lim^^oo h{s) = 0. 
Type II models: 

p 1 ( the set [—sL, sL]*^ is connected \ (t^ a\ 

[ to the boundary of [-L, Lf I ^ ^ ' ^^'^l 



L—fOo 



for any < s < 1. 



{Question: is all else ruled out? [Presumably — yes.]) 

When viewed on the macroscopic scale, models of Type I exhibit many clusters of 
still visible size, but none of them is infinite. The scaling limit of ISC in such models 
can be formulated along the lines discussed in ref. [|17| (and work in progress). 

In Type II models, the clusters visible on the macroscopic scale are qualitatively 
different than in Type I. Translated to the continuum scale, condition ( [5. 41 ) can be read 
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as saying that for any R, e > the probabihty that an e neighborhood of a given site is 
connected distance R away tends to 1, in the continuum hmit. 

The high-dimensional percolation models discussed in Section 4 are of Type II. 
One may gain a great deal of insight about their continuum limit from eq. (|4.23|) , and 
eq. ( [4 .91 ). Consider the situation in which the lattice spacing is a — >■ 0, and we observe 
the connected clusters which intersect a cubic region of fixed continuum size, r. On the 
lattice scale, the size of the box is 

W = r/a (-^ oo) . (5.5) 



Equation (|4.23|) implies that for any < 6{< d — 6) with predominant probability we 



shall see a cluster reaching distance R = a = a^^^^^^^^y^ j.id-i-S)/2 

away from the r-cube (r and R are expressed in continuum units). In fact, the number 
of such clusters diverges, according to eq. ([4.9[). If the bound produced there on Nl^w 
is a correct indication of the actual value, then not all the cluster which span the i?-box 
and intersect the r-box will intersect neighboring cubes of size r (the fraction which 
does, scales as r*^"^). The picture which is suggested by the (partial) results presented 
here is that there is a growing number of clusters reaching increasing distances, which 
in the continuum limit will typically be four dimensional, and have the graph structure 
of a tree. 

In parts, this picture is consistent with the proposal of T. Hara and G. Slade 
(mentioned in [^) that the continuum limit of large clusters looks like the process of 
Super Brownian Motion discussed in 0)- However, unlike the compact SBM, we 



expect the scaling limit to exhibit clusters of arbitrary extent. Our results suggest that 
in a certain ("weak") sense the limit, which is still to be made sense of (!), will exhibit 
two features to which we are not accustomed in lattice percolation models: 



• Percolation at the critical point. 

• Infinitely many infinite clusters. 



The latter may surprise one familiar with the general uniqueness Theorem of Bur- 
ton and Keane [^. However, it should be appreciated that the BK result requires 
discreteness (and regularity on the corresponding scale) and thus is not applicable to 
the continuum limit. 
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A. The Relation between Proliferation of ISC and Hyperscaling 



The proliferation of the Incipient Spanning Clusters is related to the breakdown 
of "hyperscaling". In order to clarify this relation, we recapitulate here the heuristic 
basis of the scaling and hyperscaling relations. Different variants of the argument can 
be found in the literature [|l], The one given below is cast in terms of quantities 

which are studied rigorously in this work. However, unlike in the rest of this paper, in 
this section we do not present rigorous results. (Rigorous results on scaling relations 
exist for d = 2 and d > d^.c. - as mentioned in Sect. 4) 



The scaling relations of some of the critical exponents can be explained by the 
following picture. The first-line assumption, related to the self-similarity, is that the 
relevant quantities scale by power laws. Next: 



1. For critical models, there are about L* clusters in with diameters of order L, 
and their volumes (defined on the U-V/lattice scale) are of the order of . If 
exceptions occur, it is assumed that they do not affect the cluster statistics by 
more than xL"^^-'. The quantities we shall look at are 

E I \'^\^ I[diam(C) > wL] j (A.l) 
\ceAi / 

with w < 1, fixed as L — oo, and k = 1,2. 

2. For p < pc there is a characteristic length 

ap) ^ {pc-pr- (A.2) 

such that: i. Only a fraction (say 1/3, or even up to (1 — L^"*^^^)) of the mean 
value of the cluster size |C(0)| is from distances greater than ^. 
a. Within the distance ^, the two point function rp(0,a;) is comparable (in the 
sense of bounded ratios) with its value at the critical point, rp^(0,x). 

Note that the mean cluster size and the two point function are in the direct relation: 
Erp(o,a;) = E,(|C(0)|) {^{pc-p)-') • (A.3) 



3. For p > Pc the density of the infinite cluster, M{p), scales as 

M{p) ^ (p-Pcf (A.4) 

and there is a characteristic length C,{p) ~ {p — Pc)~^^ such that M{p) is of the 
order of the volume-fraction of A^ which dX p = Pc belongs to clusters of diameters 
comparable with ^ (i.e.. Incipient Spanning Clusters on scale ^). 
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Let us include in the list the assumption: 



(A.5) 



which is not essential for the picture presented here of hyperscaling, but seems to incor- 
porate empirical observations, and fits well in the standard scaling ansatz. 

We shall compare now different ways of evaluating two quantities. First, let us 
look at 

5^ \C\ = J2 1[diamC(x) > wL] . (A.6) 
C c.c. in Al 
diamC > wL 



Under the assumption 3. the mean value of the expression on the right scales as L'^L^^'^+ . 
On the other hand, assuming i., the expression on the left scales as L'^L^ . Thus: 



# + D 



(A.7) 



Next, consider 
c c.c. in Az, 



\[x and y are connected ] = ^ |C(x) fl A| (A.^ 



Each of the expressions suggests a different method for evaluating the mean value. Using: 
1. for the leftmost, the definition of rj (eq. (|1.4| )) for the next, and the assumption 2. 
for the rightmost, one is led to: 



L*L 



2D 



J^'id, I jjd—2+r) jjd jjy /v— 



(A.9) 



I.e., 



# + 2D = d + 2-ri = d + -f/u.. 



(A.IO) 



A convenient rearrangement of the information (equations ( |A.7| ) and ( [A.10|) , in- 
corporating 4-) is: 



2 — 7] = / u (scaling relation) (A. 11) 

D = (dimension of the incipient spanning clusters) (A. 12) 

i^ = d- (ISC proliferation exponent) (A. 13) 
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For independent percolation the proliferation exponent # vanishes in 2D, and apparently 
also in all dimensions d < 6. In such situations we get one more equation, which exhibits 
d along the standard exponents, and hence is called "hyperscaling" : 



# = , or d- 2^ = 



(hyperscaling relation) (A. 14) 



(That brings it to four equations for six quantities.) 

Thus, the breakdown of hyperscaling is intimately related with the proliferation of 
the Incipient Spanning Clusters. However, its validity does not require uniqueness — 
just that the number of "macroscopic" clusters be typically smaller than any power of 
L, or have a finite limit in a probabilistic sense. 

Remarks: 1. It is easy to incorporate in this picture other exponents, which were not 
listed above. To make the list less incomplete, let us mention: 

Probp^ (|C(0)| > n) ^ ^ = X 

cardjccA. . diam c > ^ _ ^y-. 2-a = .(^ - #) ^^'^^^ 

(=2/5 + 7). 

where we gave the definition (for a not quite the standard one), and a scaling relation 
derived along the above lines. 

2. The above discussion is relevant also for the Ising and Potts spin models, since 
the heuristic arguments described here make equal sense in the broader context of the 
Fortuin - Kasteleyn |^8[ random-cluster models. The resulting scaling relations are the 



same in term of the recognizable characteristic exponents, except that du.c. is different 
and the significance of 7^ is lost if one is not aware of the geometric structure behind the 
spin correlations. (For Ising spin systems, hyperscaling has also another connotation: 
its breakdown implies that the scaling limit is a Gaussian random field p(]|). 
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B. Existence of Spanning Clusters in Critical Models 



In this Appendix, we prove Theorem 1, of Section 2. A key role is played by the 
following estimate. Both may be assumed to be known to experts. 

Claim: In dimension d > 1, for all s < 1/3: 

Qls - Prob, ( the set [-.L .L]^ is connected \ ^ 

^^'"^ y to the boundary of J - ^ ^ 

Proof of the Claim: By monotonicity, we may assume that 1/s is an integer. We shall 
show that if eq. ( |B.1| ) (with Cd to be specified shortly) fails for some L , then the 
connectivity function decays exponentially (at distances » L). That is well known to 
be in contradiction with the condition p = Pc TT]. 



To make the deduction, partition the lattice into cubic blocks of (linear) size sL. 
For each self-avoiding path linking a site x with y, let us associate a sequence of sL 
blocks, which are at distances approximately L apart, by the following algorithm. The 
zero-th block is the one containing x. Next is the block which the path reaches when it 
hits the boundary of the cube of size (4 + s)L concentric with the first block, and so on: 
once a stopping point and a block are selected, we center on the sL block a large cube 
of size (4 + s)L, and let the next stopping point be the exist site, and the next block be 
the corresponding element of the lattice cubic partition. Notice that for each block in 
this sequence, other than the end points, the event seen in eq. ( |B.1| ) occurs twice (at the 
end points once) and disjointly so. The van den Berg - Kesten inequality |31|, which is 



described above, permits to deduce that the probability that x and y are connected by 
a self-avoiding path which corresponds to given sequence of {k + 1) blocks, is < Q 



2K 
L,s' 



Summing over the possibilities we get: 

t{x, y) = Probp^ ( x and y are connected 

< Const. e-^\^'y\ , (B.2) 



< E 

k>\x~y\l{2L) 



with /X > if 2d 5^"^ 4^ Q\ ^ < \. Since /x = at the critical point (|5l|, |ll| 
deduce that eq. (E^) holds, with C^ = V2d 5'-'^~^^/^. 



Proof of Theorem 1: If the probability that that the slab S^^t is traversed is very 
small, then so will be Ql/2,s ^ot s = 1/2 — t. To quantify that, let us consider two 
concentric cubes, [-L/2,L/2]'^ and [-sL/2, sL/2]'^. If none of the 2d similar slabs which 
enclose the smaller cube is spanned (in the short direction) then that inner cube is 
not connected to the outer's boundary. Since the 2c? events are positively correlated 
(IPPj, 1^), the probability of their simultaneous occurrence is greater than the product. 



Thus: 

> l-RL,t. (B.3) 
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The combination of eq. ( |B.3| ) with eq. ( [B.l| ) shows that for < t < 1/2 the spanning 
probabihty is positive as claimed in eq. with h{t) > Cd (1/2 — t)^'^^^^^'^. 



The crossing of narrow slabs, with t \ 0, can be estimated by cutting the slab 
into {3/tY~^ smaller ones with the aspect ratio close to 1/3. Their spanning proba- 
bility is uniformly positive ii p > Pc (as shown by the previous discussion). Using the 
independence of the events, one obtains eq. (|2.2|) . 



C. Uniqueness for Supercritical Models 

In this appendix we supplement Theorem 2, by proving that only for p = Pc would 
there be positive probability for observing more than one spanning cluster in arbitrarily 
large systems. We shall use the proven fact that in d > 2 dimensions pc is the limit of the 
quadrant-percolation thresholds for slabs of finite width |53|. This is an important 
technical statement, which makes a variety of 2D arguments applicable to supercritical 



models in dimensions d > 2 ||2^. (In particular, it permits to prove that the spanning 



probability itself tends to 1, as L ^ oo.) 



Theorem 6 For any p ^ p^, and t > 0, 

I 

DL{t,p) = Probp 



there is more than one spanning 
y cluster m SL,t{ = [0,tL] x [-L, L]('^-i)) 



(C.i: 

L^oo 



Proof: For any p < p^ the two point function decays exponentially |5^, and that 



easily implies: DL{t,p) < Const. L — >■ (for L — >■ oo). The assertion which still 
requires a proof is that for p > Pc (and d > 2) the probability of there being more than 
one spanning cluster tends to 0. 

Let us define, for pairs of sites on the left boundary ({xi = 0}) of the semi-infinite 
cylinder [0, oo) x [0,1]'^'^: 

„ , ^ 1 ( X and y are in distinct spanning clusters \ 

^^^^'^) = ^^^^4 of[0,4x[0,L]-i )■ (^-2) 

The events seen in eq. ( |C.2| ) are obviously decreasing in fc, and thus the ratio 
Gk+w{x,y)/Gk{x,y) represents a conditional probability. For p > pc, let w{p) be the 
smallest slab width for which there is percolation in the quadrants [0,w]'^^^ x [0, oo)^ 



w{p) < oo p2| , [53[|). The above conditional probability, of the {k + W^)-th event con- 



ditioned on the k-th. event, is uniformly smaller that one, since there is a uniformly 
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positive probability that the two distinct clusters will be joined by a path within the 
added slab. (The reason is explained more explicitly in Section 4 of ref. in the 

context of a rather similar argument.) Thus 

Gk{x, y) < Gk-^{x, y)e-'' ... <A e-^^'^ (C.3) 

with some a > 0. Consequently: 

_ , / there is more than one \ ^ ^ / n 

^ > spannmg cluster m ' ^ 



x, y e d[0,2LY 

xi, yi = 

< A (2L)2'^e-"2^/"' — > . 



(C.4) 



Remark: We see that the probability of there being more than one spanning cluster is 
exponentially small, in L, for p both above and below Pc, but not for p = Pc (eq. ( p.3|) ). 
This quantity yields a natural and meaningful characteristic length scale ^{p), and could 
perhaps offer a useful tool for a more thorough analysis of the the scaling relations. 
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